Geometrical spin dephasing in quantum dots 
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We study spin-orbit mediated relaxation and dephasing of electron spins in quantum dots. We 
show that higher order contributions provide a relaxation mechanism that dominates for low mag- 
netic fields and is of geometrical origin. In the low-field limit relaxation is dominated by coupling 
to electron-hole excitations and possibly 1// noise rather than phonons. 
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Recent experiments 0, 0] demonstrate that spins of 
single electrons confined in quantum dot structures can 
be manipulated in a quantum coherent way, thus opening 
exciting perspectives for quantum information processing 
Thus a thorough understanding of spin relaxation 
and decoherence processes is crucial. This requires identi- 
fying the sources of fluctuations and the mechanisms how 
they couple to the spins, as well as analyzing the non- 
equilibrium decay laws. So far, two major mechanisms 
of dephasing have been identified: One is the hyperfine 
coupling to randomly oriented nuclear spins While 
in a free induction decay this leads to fast dephasing, 
on a time scale of order ~ 10 ns, the effect of the quasi- 
static nuclear field can be eliminated largely by spin-echo 
techniques, bringing the dephasing time into the range 
of ~ Ifis ',2|. The second mechanism is the coupling to 
(piezoelectric) phonons in the presence of spin- orbit in- 
teraction Phonons create a fluctuating electric 
field acting on the electrons' orbital degrees of freedom, 
which couples via the spin-orbit interaction to the spin. 
If time reversal symmetry is broken by a magnetic field 
B the process leads to spin relaxation. This mechanism 
is important in sufficiently strong fields 5], however, as 
usually described in the literature, it is ineffective in van- 
ishing fields for electrons in the ground state doublet. 

In this paper we study, what processes destroy the spin 
coherence in vanishing or low magnetic fields. We show 
that higher order (e.g., two-phonon) virtual processes, 
usually neglected in the literature, provide a relaxation 
mechanism that persists as B — *■ 0. These relaxation 
processes are of geometrical origin and related to the dif- 
fusion of the Berry phase. Berry phase emerging from the 
spin-orbit interaction was introduced in Ref . ^] . A differ- 
ent phenomenon also called geometric dephasing was dis- 
cussed in Ref. . The processes we consider are the ana- 
logues of Elliott's spin relaxation in bulk semiconductors 
and metals [T3 |. for which a geo metrical interpretation 
has recently been given in Ref. [Ill • Processes of this type 
have also been studied in the context of phonon-induced 
relaxation in electron spin resonance experiments |l2l |. 
and were found to lead to a relaxation even as i? ^ 0. 



However, the geometric origin of the mechanism has not 
been revealed. Moreover, the truncation of the Hilbert 
space to two orbital levels, used there, is insufficient, since 
amplitude cancelations from higher orbitals prove to be 
crucial. 

Moreover, we observe that spin relaxation is induced 
by any kind of electric field fluctuations, not merely by 
phonons. In fact, in order to confine, control, and mea- 
sure the electron one attaches electrodes and quantum 
point contacts to the qubit They produce Ohmic 
fiuctuations with dominant spectral weight at low fre- 
quencies. For typical quantum dots we find that such 
Ohmic fluctuations provide the dominant mechanism for 
Berry-phase dephasing at i? = 0, and they also provide 
the leading channel for spin relaxation at fields below 
roughly 1 Tesla. In addition to Ohmic fluctuations, the 
quantum point contacts produce shot noise when driven 
out of equilibrium, which further relaxes the spin |l3l |. 
Finally, 1// background charge fluctuations, present in 
most mesoscopic systems, also couple to the orbital mo- 
tion of the electrons and dephase the spin. In the follow- 
ing we shall develop a formalism that allows us to treat 
different types of environments on an equal footing, and 
also to take into account higher order virtual processes, 
leading to a geometrical dephasing. 

We consider a single electron confined to a lateral 
quantum dot by the potential V^(r), in the presence of 
a magnetic field B. To be specific, we assume the field to 
be oriented parallel to the plane of the dot, but our pro- 
cedure can be generalized to arbitrary directions. The 
static part of the Hamiltonian then reads 



Hn = 



P^ + Vir)-'J^B-a 



HsO = Ol[pyax - Px(7y) + P[Py<yy - Px<yx) ■ (2) 

The magnetic field couples to the electron only through a 
Zeeman term with ^-factor g |Q| . The last terms describe 
the Dresselhaus (/?) and Rashba couplings (a) between 
the spin a of the electron and its momentum 0, Il5l |. 
For a dot of size d the typical energy of the spin-orbit 
coupling scales as ^ /3/d, while the level spacing scales 
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Here 5Vint (i) denotes the fluctuating part of the Hamil- 
tonian in the interaction representation and P = 
X](T We separated terms that involve direct tran- 

sitions between the two lowest states from transitions 
via excited states. In the spirit of an adiabatic approxi- 
mation, these latter processes can be integrated out to 
yield an effective Hamiltonian in the two-dimensional 
subspace. Technically, this is performed by introducing 
slow and fast variables, t = {ti + t2)/2 and t = ti — t2, 
in the last term of Eq. Q), 



FIG. 1: Left: Single electron lateral quantum dot in a mag- 
netic field, which lifts the ground state degeneracy. Virtual 
transitions to excited states are induced by weak fluctuations 
of the external fields 5Ex{t), 6Ey{t). Right: Graphical rep- 
resentation of the evolution operator. Virtual transitions to 
excited states n (wavy lines) are integrated out to yield an 
effective Hamiltonian within the doublet subspace. 



as ujQ ~ l/{m*cP). Therefore, for dots with small level 
spacing, uq < 1 K, the spin-orbit coupling cannot be 
treated perturbatively. 

Next we account for time-dependent fluctuations of the 
electromagnetic fleld, which add a term 



-it(e„-£„/)-ir[i(e„+e^,)- 



SV^n{tl)SVna'{t2) 



SV{t) = X^{t) 6^ 



(3) 



to the Hamiltonian. (A summation over repeated in- 
dexes, such as /i, is assumed throughout.) The terms 
denote independent operators in the Hilbert space of 
the conflned electron (e.g., x, y, x^, . . .), while the terms 
X^(t) denote the corresponding fluctuating (in general 
quantum) fields (e.g., SE^, 5Ey, Vx^Ex, ...). They may 
be generated by various environments, such as phonons, 
localized defects, or electron-hole excitations. Informa- 
tion about their specific properties is contained in the 
spectral functions, to be specified later. Note that this 
formulation covers also quadrupolar fluctuations. 

In case of time-reversal symmetry the ground state of 
the dot is two-fold degenerate. This degeneracy is split 
in an external magnetic fleld. If g^sB <^ luq, and as 
long as the noise is adiabatic with respect to the orbital 
level splitting, T <^ luq, the dynamics of the spin remains 
constrained to these two states. Under these conditions 
(following the method described in Ref. Hii) we can derive 
an effective Hamiltonian for the two lowest eigenstates 
\a = ±) by expanding the evolution operator U{t) — 
T exp{— i Jg dt' 5Vint{t')} and projecting to the subspace 
{|(t)} as sketched in Fig. ^ This yields 

PU{t)P = 1 - « / dtiPSVint{tl)P 

Jo 

- f f PSVi,,t{ti)PSVi,,t{h)P (4) 

J Jti>t2 

PSVintiti)il^P)SV,ntit2)P + -- - ■ 

Iti>t2 



expanding the interaction potential in r as 5V{ti 2) ~ 
SV{t) ± ^ ^'^^(0 + • ■ -I and integrating with respect 
to r. Here Eo- and e„ denote the eigenenergies of the 
lowest doublet and higher eigenstates of Hq, respectively. 
In this way the last term in Eq. Q becomes local in 
time. Retaining only processes up to 2nd order, we find 
an effective Hamiltonian within the lowest-energy two- 
dimensional subspace, characterized by the 'pseudospin' 
Pauli matrices Tx,y^z, 

+ ^(^X^^X" -X^^X^^Cj^J -f. (5) 

Due to the spin-orbit coupling, which is not assumed to 
be weak, eigenstates do not factorize into orbital and 
spin sectors (hence the term 'pseudospin'). The static 
effective field, Bctt = (e+ — £-)z, accounts for the spin- 
orbit renormalization of the g-factor and defines the z 
direction in the doublet space. The couplings C^'^\ de- 
termining the effective fluctuating magnetic fields felt by 
the pseudospin, are given by 
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(6) 
(7) 

(8) 



The summation is restricted to excited states of higher 
doublets n ^ a, a'. We do not provide exphcit ex- 
pressions for the eigenenergies e^, £n, matrix elements 
O^^/, Offn, or couplings C^^\ but below we will evaluate 
them numerically and provide quantitative estimates for 
a generic model. We further note that both C?^^"* and 
Cij,J turn out to be transversal to Bcs, therefore con- 

-^(2) 

tributing only to relaxation, whereas C^iJ has in general 
also a parallel component that leads to pure dephasing. 

In time- reversal symmetric situation, (i.e. for B = 0), 
the first three terms of Eq. ||SJ| vanish identically Q . Only 
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the last term survives, and leads to spin dephasing. It 
has a geometrical origin. To demonstrate this, let us as- 
sume that the fluctuating (adiabatic) fields are clas- 
sical. We introduce the instantaneous ground states of 
the Hamiltonian, = \^n{Xf^{t))) defined through 

the equation 



[Ho + SV{X^)]\<S>^{X^,)) = E^{X^)\^^{X^)) 



(9) 



Noting that, to lowest order perturbation theory, the two 
degenerate instantaneous ground states are simply given 

by « \a) + j:n'HH^VW)/i^- - £") ' can 

rewrite the last term in Eq. jSJ in the familiar form 



HfAB ^ 0) ^ -^{^\^., 



(10) 



which shows clearly that the last term is due to a gen- 
eralized (possibly non-Abelian) Berry phase [tH. fl3 Il9l| 
acquired in a degenerate 2D subspace. In vanishing mag- 
netic field, Eq. H10() can be shown to hold to all orders of 
perturbation theory within the adiabatic approximation. 
If at least two linearly independent fluctuating fields cou- 
ple to the dot, they can produce a random Berry phase 
for the system and cause geometric dephasing at i? = 0. 
When more noise components are present, the Berry 
phase may become non-Abelian and all components of 
the spin may decay. 

So far, our treatment has been rather general, ap- 
plicable for arbitrary noise properties and dot geome- 
tries. In its full glory, Eq. ||SJ) describes the motion of 
the pseudo-spin coupled to three fiuctuating "magnetic 
fields" . In general, the dynamics induced by these non- 
commuting fields is complicated. To obtain a qualitative 
understanding of the dynamics we analyze the spin re- 
laxation and pure dephasing times |2Cl| . Ti and Tj* (with 
I/T2 = l/2Ti + 1/T^). They are defined only for suf- 
ficiently strong effective fields, BcS > 1/T2*- 
the limit B = we evaluate what we call the geometrical 
dephasing time Tgoom- For the quantitative estimate we 

consider a parabolic confining potential, V{r) = |rp 
with level spacing luq and typical size xq — 1/ ^ujQun* . 
Furthermore, we take into account only dipolar fiuctu- 
ations, X = eSE^xa and Y = eSEyXo coupling to the 
operators Ox = x/xq and Oy = v/xq, respectively. 
We assume the two components X and Y to be inde- 
pendent of each other, but with identical noise spectra, 
Sx{i^) = Syito) = S{u!) = 7rg(a;)coth(w/2T), with g^uj) 
being the spectral function of the bosonic environment 
(phonons or photons). 

The spectral function g{(^) for phonons can be esti- 
mated along the lines of Ref . |^ For the parameters spec- 
ified in Ref. we find for piezoelectric phonons in typi- 
cal GaAs heterostructures at low frequencies, f»ph(a-') = 
xl Aph uj^ with Aph « 4 • lO^'^K-^nm-^ 14J. With 
these parameters we obtain relaxation rates generated 
by the first term in Eq. jSJ that coincide with those of 



Ref. at not too high values of the field [ij. Similar 
values are obtained for the parameters of Ref. For 
Ohmic fluctuations the spectral function is linear at low 
frequencies, gnioj) = An w jU 0|. The prefactor An 
depends on the dimensionless impedance of the circuit. 
An X R-^I^]- For typical values of the sheet resistance 
of the 2-DEG (10^ - lO^Q/a) we estimate it to be in the 
range 0.1 > An > 0.01. For 1/f noise the power spec- 
trum is S{lu) ~ Xi^f/\uj\. We will further comment on 
its strength below. 

We first estimate the contributions T^^^ and T2**'*\ 
derived from the three terms (i = 1,2,3) in Eq. jSJ, 
for a non-vanishing in-plane magnetic field, B^g 3> 
l/Ti, I/T2*. The coupling C^^^ turns out to be per- 
pendicular to BcS 0, and for low magnetic fields and 
weak spin-orbit coupling is proportional to jC^^^I ~ 
max{a,/?}. This fluctuating field therefore con- 



tributes to the Ti-relaxation only. 
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(1) 



= 2 \C 



Sx{Bcs) 



(11) 



It scales as l/T^^^-* ~ B^ max{i3, T} for Ohmic dissipation 
and as ^ B"^ m.a.x{B,T} for phonons. As a consequence, 
for dots with level spacing in the range uiq w 1 ... 10 K 
Ohmic fluctuations dominate over phonons for low fields 
with B < 1 . . . 3 T. 

The second term, IC^^^I ^ max{a^,/9^} B gives 
rise to both relaxation and dephasing. The two rates are 
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(2) 



- 4 + + 2|4^^ip) Sxy{B.^) , 

' =4(|c(,^i?J^ + l4^vlSP + 2|4^?J^)^xW0), 



■(2) 



/„/ UJ-\-LJ \ U)—U) \ 
1 -cosh(c;;/2T)/cosh(w/2r) ' 

with C'/ii',i"^"'' denoting the symmetrized component of 
Cfiil perpendicular/parallel to -Boff ■ Thus, for Ohmic dis- 
sipation l/Tp' vanishes as ~ maxji?"^, T'^}, while for 
phonons it scales as ^ B^ max{i?^, T^}. 

C'^^^ is also perpendicular to Beg. Its contribution to 
the relaxation is 



1 



(3) 



(12) 



S 



XY-XY 



{to) 



~,2 „( U>+LU \ „f u — il] \ 



1 - cosh(c;)/2T)/cosh(cj/2r) ' 



with CfiiJ,a being the anti-symmetrized compo- 
nent of ■ Most importantly, with ICj^yl ~ 
■i 4 maxjo;^, the rate l/T,'^' approaches a non- 
zero value at low fields, \/Ti,\/T:^ <^ B <^ T, and 
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FIG. 2: Spin relaxation rates for a GaAs quantum dot 
with level spacing uuo = IK as function of the Zeeman field. 
We chose a temperature T — 100 mK and Ohmic coupling 
An — 0.05. Below B* « 1 T spin relaxation is dominated 
by coupling to Ohmic fluctuations. For B < B** « 15 mT 
geometrical spin relaxation due to coupling to Ohmic fiuctua- 
tions dominates. For all B values plotted, the Bloch-Redfield 
consistency requirement, BeS 2> 1/T^^\ is satisfied. 

scales as ~ niax{_B^,r^} for Ohmic dissipation and 
~ max{_B^, T^} for phonons. 

Finally, at i? = the geometric dephasing rate is given 
by 1/T^^\ extrapolated to zero field 

= 2\CPyJ'S^y_^y{u; ^ 0) . (13) 

gcom 

In our example with only two noise components this pro- 
cess dephases only the components of the spin perpen- 
dicular to C^y. 

The relaxation rates corresponding to the different 
terms in Eq. ^ and various noise sources are shown in 
Fig. H Clearly, for external fields B < B* « 1 . . . 3 T 
Ohmic fluctuations provide the leading relaxation mech- 
anism. The crossover fleld B* is not very sensitive to 
the specific value of the spectral parameter An and is 
independent of the spin-orbit coupling. Below a sec- 
ond crossover fleld, B** w 15mT, the geometric de- 
phasing induced by Ohmic fluctuations starts to dom- 
inate. This second crossover scale is very sensitive 
to the spin-orbit coupling and temperature, scaling as 
B** ^ max{Q!,/3}(l/a;o)(T/a;o)^. E.g. for a level spacing 
ujo - IK and temperature T = 100 mK (T = 50 mK) 
the Berry phase mechanism gives a relaxation time of 
the order of 700 /is (20 ms). For even lower tempera- 
tures or smaller dots with level spacing cuo ^ lOK the 
-B — > relaxation time is quickly pushed up to the range 
of seconds. 

Finally, we comment on the effect of 1// noise. In 
most cases, the non-symmetrized correlators for 1/f 
noise, needed to calculate correlators as Sxy or Sxx, 
are not known. Yet, for \lu\ <^ T we can provide an 



T 

estimate Sxy{^^) ~ / §^Sx{oJ — Ct')S'y(w), and sim- 

-T 

ilarly for S^y-XY- '^^'^ B — geometrical dephas- 
ing rate due to the 1/f noise can be estimated as 
T^clm ~ l<5xFpA^/^(r) Wc, where tUc is the upper fre- 
quency cut-off for the 1// noise . Accounting for the 
high-frequency (Ohmic) noise, sometimes observed to be 
associated with the 1/f noise HI HI El 113 , the esti- 
mate becomes Tg-^^ « |C^^|-|2a2/^(T) T. While the 1/f 
noise of background charge fluctuations is well studied 
in mesoscopic systems, the amplitude of the 1/f noise 
of the electric fleld in quantum dot systems is yet to 
be determined. If we assume that this noise is due to 
two-level systems at the interfaces of the top gate elec- 
trodes, we conclude that in the parameter range explored 
here, the effect of 1// noise is less important than that 
of Ohmic fluctuations. However, in quantum dots with 
large level spacings in the low-temperature and low-fleld 
regime, these fluctuations could dominate over the effect 
of Ohmic fluctuations and eventually determine the spin 
relaxation time. 
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